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Abstract. The idea that geomagnetic westward drift indicates convective leveling of the planetary momentum 
gradient within Earth’s core is pursued in search of a differentially rotating mean state, upon which various 
oscillations and secular effects might be superimposed. The desired state conforms to roughly spherical boundary 
conditions, minimizes dissipative interference with convective cooling in the bulk of the core, yet may aid core 
cooling by depositing heat in the uppermost core and lower mantle. The variational calculus of stationary 
dissipation applied to a spherical vortex within the core yields an interesting differential rotation profile, akin to 
spherical Couette flow bounded by thin Hartmann layers. Four boundary conditions are required. To concentrate 
shear induced dissipation near the core-mantle boundary, these are taken to be: (i) no-slip at the core-mantle 
interface; (ii) geomagnetically estimated bulk westward flow at the base of the core-mantle boundary layer; (iii) no- 
slip at the inner-outer core interface; and, to describe magnetic locking of the inner core to the deep outer core, (iv) 
hydrodynamically stress-free at the inner-outer core boundary. By boldly assuming the axial core angular 

momentum anomaly to be zero, the super-rotation of the inner core relative to the mantle is calculated to be at most 
1.5°/yr. 

1. Introduction 

The westward drift of certain features in the geomagnetic field across Earth’s surface at a few tenths of a 
degree per year has long been used to suggest westward drift of Earth’s magnetized core (Halley, 1692). To the 
extent that the field is frozen into an electrically conducting fluid outer core, westward drift indicates westward flow 
near the top of the core. Bulk westward flow at about 0. 12°/yr emerges from somewhat more detailed inversions of 
secular geomagnetic change (see, e.g., Voorhies, 1995); however, this rate can be reduced by diffusion, hence 
slipping, of the field through the imperfectly conducting fluid (see, e.g., Voorhies, 1993). Such westward flow poses 
an interesting problem for those who suspect that, during the course of geologic time. Earth’s core and mantle have 
relaxed to a state wherein the angular momentum of the core is quite close to that of a rigid body rotating with the 



angular velocity of the mantle. (This apart from differential spin up due to inner core solidification, differential tidal 
spin down, differential precession, and various oscillations). 

A slight convective leveling of the planetary azimuthal momentum gradient still seems to offer the most facile 
explanation of westward drift. Compared with a non-convecting, uniformly rotating state, buoyant, low azimuthal 
momentum fluid at depth floats towards the upper surface and drifts west; dense, high azimuthal momentum fluid 
above sinks towards the lower surface and drifts east; and the azimuthal momentum gradient is reduced. The 
reduction is arguably very small for several reasons: the convective motions are feeble, strongly constrained by the 
rotation itself, and may be largely confined to spirals on surfaces of constant planetary momentum. Slight 
departures from motion on such surfaces may nonetheless be realized, notably if they increase convective efficiency 
by decreasing dissipation or, perhaps for thermo-compositional magneto-convection, by distributing Ohmic and 
viscous dissipation so as to aid core cooling. Westward drift has thus long provided a reason to consider seriously 
the convective outer core geodynamo hypothesis, whereby convective fluid motions maintain the aperiodically 
reversing geomagnetic main field against magnetic diffusion and Ohmic dissipation. 

This facile explanation of westward drift clearly requires eastward flow of the fluid at depth in the outer core. 
If the total angular momentum of the core is quite close to that of a rigid body of identical inertial moment rotating 
with the mantle, then the speed of the eastward flow at depth will exceed that of the westward flow above so long as 

the geometric effects of reduced moment arm and volume at depth dominate the increased density of quasi- 
hydrostatic compression. 

The foregoing argument for eastward flow at depth has long been obvious to students of the secular variation. 

It is now clear that tight magnetic coupling of eastward flow near the base of the outer core to the similarly high 
conductivity solid inner core would tend to spin up the inner core. Indeed, deviation from co-rotation shears the 
magnetic field lines that thread the inner core boundary and are thus embedded in both liquid and solid conductors. 
Such shear induces electric current and restoring Lorentz torques that tend to return the system towards a 
magnetically locked equilibrium - perhaps after several resistively damped oscillations (Gubbins, 1981). The 
essential geophysics was elucidated by Glatzmaier & Roberts (1995a, b; 1996), via numerical simulation of the core 
geodynamo, and led to their celebrated prediction of an eastwardly drifting inner core. 

Seismologic substantiation of this prediction by Song & Richards (1996) and Su et al. (1996) suggests 
considerable uncertainty as to the amplitude of eastward inner core drift. The former obtain about l.l°/yr; the latter 



obtain about 37yr (2.27+0.90 °/yr, 3.02 + 0.427yr, and 3.287yr depending upon details of the inversion). Creager 
(1997) obtains a lower rate of about 0.257yr, while Souriau et al. (1997; 1998) question the seismologic detection of 
inner core rotation itself. This paper shows how eastward drift of a solid inner core can be simply, albeit not 

unambiguously, calculated from geomagnetic westward drift without recourse to either numerical simulation or 
models of inner core acoustic anisotropy. 

2. Simple Model and Example 

Consider a planet rotating with sidereal angular velocity ft(r,t) and mass density equal to marginally stable 
reference density p(r) plus small perturbation 8p(r,t) ( 1 8p/p | « 1). For simplicity, position r is measured in 
centered spherical polar coordinates (radius r, colatitude 0, and east longitude <j>) rotating with a rigid mantle at 
uniform angular velocity Cl Q z (z = rcosG - §rsin0). This filters out secular tidal despin, polar wander, precession, 
nutation, and decadal and higher frequency fluctutations in the angular velocity of the mantling solid, which fall 
outside the focus of this paper. 

The planetary rotational velocity £> 0 zxr is eastward Q Q rsin0$. The planetary momentum density M i s 
eastward M^, is 4>p£2 o rsin0, and is the planetary azimuthal momentum of a material parcel of unit volume. Clearly 
tends to increase with distance away from the rotation axis, s = rsin0. For a small planetoid of homogeneous 

mass density p 0 , the planetary momentum gradient VM^ is uniformly steep p o n Q s. For a self-gravitating 
compressible planet, 

VM(j> = pQ 0 (s + sVp/p) . (1) 

If the rotation is not too fast, then a roughly spheroidal planet remains nearly spherical, the stratification is nearly 
radial p(r), Vp ~ rd r p, and 3 r p < 0. If the density scale height is not too small (p I d r p i > s), then VM^ is not 
vastly different than for homogeneous density; however, surfaces of constant planetary momentum flare away from 
the rotation axis as distance from the equatorial plane I z I increases. Such surfaces are more akin to hyperboloids of 
revolution than right circular cylinders. For terrestrial example, calculated from the seismologic density model of 
Kennet et al. (1995), the constant surface tangent to the equator of Earth’s solid inner core intercepts the fluid 
outer core-mantle boundary at colatitude 25.4° - some 300 km south of the right circular tangent cylinder intercept 


at 20.5°. 



To avoid convective leveling of the planetary momentum gradient in a thick, rapidly rotating, roughly 
spherical annulus such as Earth’s electrically conducting liquid outer core, one need only arrange (a) perfect 
exchange of planetary momentum between rising and falling fluid; (b) perfect confinement of convective motions to 
surfaces of constant planetary momentum; or (c) perfect exchange of planetary momentum between fluid parcels 
moving onto and parcels moving away from each surface of constant planetary momentum. Although case (c) is 

thought to be an excellent approximation, it seems unlikely that any of these conditions could be perfectly satisfied 
in a real planet. 

Core Angular Momentum Anomaly. The angular momentum density in the core is the sum of planetary and 
deviatoric portions. The planetary angular momentum density is northward (rx^M^ = -§pf2 Q r 2 sin0). The 
deviatoric angular momentum density is due to perturbation density and relative velocity v(r,t). The total angular 
momentum of the core is the volume integral 

L = CV Cr4pa ° rSinG + rX(pV + 8pv + ^o rsin0 )] dV (2a) 

s Lp + AL , 

where CV denotes core volume. If p(r) is axisymmetric p(r,0), then the planetary angular momentum of the core L 

P 

is parallel to the reference rotation axis. The deviatoric angular momentum, 


AC = / J J rx(pv + 5pv + <J>8pfLj-sin0) dV 
CV 


~ J J J rx(pv + <)>5pQ rsin0) dV , 
CV 


is just the core angular momentum anomaly which some think small ( I AL I « I L p I ). 

Does the Taylor-Proudman Theorem Inhibit Leveling? To find out, partition the mass transport equation 
into main (V*pv = 0) and fluctuating parts. The curl of an approximate balance between the Coriolis force density 
(2Q 0 zxpv) and hydrostatically uncompensated scaloidal force densities (such as perturbation pressure) gives the 
Proudman-Taylor constraint 2Ct 0 (£.V)pv « 0 (see, e.g., Gubbins & Roberts, 1987). The relative momentum density 
pv does not vary much along the direction of the rotation axis; yet motions orthogonal to the rotation axis are 
allowed and might slightly level VM^. Such leveling would appear as a non-zero mean gradient in the deviatoric 
azimuthal momentum (pv^ + Spv^ + 5pQ Q rsin0) rather than a change in VM^ itself. 



Purely azimuthal flow v^Cr) satisfies the constraint when relative momentum pv^ varies only with s (i.e., 
V r) ” fP( r )l _1 [f(rsin0) + K 0 ] where K Q is a constant and f is a function of s alone). At fixed s, the decrease in p 
with radius, hence distance I z I from the equatorial plane, can be compensated by faster flow. Surfaces of constant 
> thus tend to curve towards the rotation axis as f z| increases. Such curvature resembles a spherical vortex 
modulating azimuthal flow that would otherwise be invariant along more familiar co-axial right circular cylinders. 
The effect could be quite large for motions spanning multiple density scale heights, as in a stellar or giant planetary 
convection zone, but the mass density contrast across Earth’s outer core is only 22%. The notion of a spherical 
vortex persists because Taylor-Proudman conditions are not fully met by magneto-convection of a viscous fluid; 
because smaller scale motions might excite and maintain such a vortex; because of the shape of the bounding solids; 
and because of geomagnetic westward drift. 

2.1 Simple Spherical Vortex and Boundary Conditions 

As in many models, inner-outer core and core-mantle transition zones are here approximated by sharp material 
interfaces. The solid inner core is approximated by a rigid ball of radius d, the fluid outer core by a spherical 
annulus of outer radius c and inner radius d, and the reference density by laterally homogeneous p(r). For the bulk 

eastward or westward flows considered here, spherical shells rotate as if rigid bodies; therefore, attention is directed 
towards the spherical vortex 

V(r) = w o (r)sin0$ = o)(r) rsinB $ , 

where w Q (r) is the equatorial amplitude of bulk eastward flow, to(r) is the perturbation angular velocity or "eastward 
drift", and the relative vorticity is 2cn(r)S. Spherical vortex (3a) is partitioned into frozen field and diffusive portions 

co(r) = tOjjfi) + ayr), (3b) 

where to B represents geomagnetic drift in the frozen-flux approximation and to^(r) denotes diffusive slip. The 
radially varying ratio I to/to^ I indicates the relative importance of advection to diffusion in the vortex. 

The hypothetical kinematic model (3a) is not intended to satisfy conditions needed for the Taylor-Proudman 
theorem, the rigidized sphere model of differential precession (Vanyo, 1984), a dynamical simulation (Glatzmaier & 
Roberts, 1995a, b; 1996), or a thermal wind model (Aumou et al„ 1996; 1998). The rigidized sphere model would 
correspond to a homogeneous spherical vortex beneath a boundary layer (tn(r < c') = 0.1°/yr). Whether or not an 
inhomogeneous spherical vortex is excited and maintained by smaller scale, possibly turbulent, convective and/or 



wave motions (u - v - V), it conforms to the boundaries more naturally than a cylindrical vortex and thus eases 
discussion of boundary conditions. 


Although the outer core might slip past the mantle almost as easily as mercury under leaded glass, radius o' is 

used to denote the base of a thin viscous boundary layer separating a free-streaming fluid core from a rigid mantle. 
At the edge of the free-stream, 


®(C) = cOgCc') + ©^(c') 

0 ) B (c‘) = -0.11577yr. <3c) 

The illustrative numerical value is the 1945-1980 mean from the preferred solution of Voorhies (1995; the value 

-0.1217yr was for 1967.5). Inversions for steady flow and steady flux diffusion suggest some eastward diffusive 

slrp 0.5 1 0) B (c-)l > eye) > 0; steady, surfically geostrophic flow inversions suggest a westward diffusive slip 

“n (C ) £ 0 and slower w «stward magnetic drift (Voorhies, 1993). The maximum westward flow used below is the 
frozen-flux extreme, -co(C) = -cog(c"). 

For a rigid mantle, the hydrodynamic no-slip core-mantle interface condition is 


w(r>c) = 0. 

For a rigid inner core, the hydrodynamic no-slip inner core boundary condition is 


(3d) 


co(r < d) = ©(d) = lim ©(d + 5) . 

8 -> 0 + (3e) 

Viscous and magnetic locking of a rigid inner core to fluid at and near the base of the outer core suggest little or no 
slip, and little or no shear, at the inner boundary. At the desired equilibrium, neither advection nor shear of 
magnetic field lines threading the inner core excite restoring Lorentz torques. The no shear condition amounts to the 
hydrodynamic stress-free inner boundary condition, 

co'(d) = 0, 

(3f) 

where © denotes d r w. The combined no-slip/no-stress conditions yield ©(d‘) = ©(d) = w(d + ). The intent of (3f) is 
to filter out, not rule-out, oscillations about, and secular trends in, the desired equilibrium. Granting conditions (3c- 
0, diffusively uncompensated westward drift to(c-) - co B (c') offers kinematic driving for a non-trivial spherical 
vortex. Gravitational locking of the inner core to the mantle (Buffet, 1997) would further require to(d) to be zero. 



2.2 Axial Angular Momentum Anomaly for Small Density Perturbations 

For spherical vortex (3a), the integrand in equation (2c) for AL Z would vanish everywhere if Sp/p were equal 
to - 0 >(r)/n o . In the frozen-flux limit of (3c), such detailed balance would have 8p(c')/p(c‘) be 8.8xI0‘ 7 . Provided 
the volume integral of Sp vanishes, this density surplus above would imply a density deficit, hence eastward flow, 
below. For example, if suiplus mass in the shell of radius c‘ and thickness Sr were balanced by a mass deficit in an 
equally thin shell of radius d + , then Sp(d + ) would be -(c/d) 2 Sp(c-), or -8.2 6p(c ); moreover, with e * d/c, 

co( d+ ) = -n o 8p(d + )/p(d + ) = I2 0 [c/d]" 2 8p(c‘)/p(d + ) = -<o(c‘)e 2 p(c‘)/p(d + ) ( 4 ) 

would be -6.67m(c-) (using e = 0.3499 and p 0 (O/p 0 (d + ) = 0.8168 from Kennetet ah, 1995). The cotresponding 
eastward drift of a magnetically locked inner core would be about 0.8°/yr. 

If buoyancy and Conolis force densities are the same order of magnitude in a convecting outer core, then 5p/p 

" l 2 £ 2 0 xv/gl = 4x 10 ' 9 . The density perturbation needed for detailed balance, being up to 200 times this value, 

seems quite unstable. Perturbation angular momentum density shall thus be omitted compared with relative angular 

momentum density (rxpv). The axial angular momentum anomaly (2c) for spherical vortex (3a) thus reduces to 
2k n c 

AL Z = J J J rsinG [p(r)w(r)sin0] Ain0 drd0d<b ^ 

0 0 0 (5a ' 

c 

= (8n/3)/ p(r)oXr)r 4 dr . ^ 

Although p(r) has long been inferred seismologically with great confidence, o>(r) is needed to evaluate AL 
From (3c), the key supposition AL Z = 0, and several arbitrary differential rotation profiles and boundary conditions, I 
calculated inner core drifts ranging from 0.2°/yr to 37yr. The agreement with the range of seismologic estimates 
appears wholly fortuitous. The following sample reveals the ambiguities and shortcomings of such ad hoc 
calculations, yet suggests a more satisfactory approach. 

2.3 A Parabolic Shot at Inner Core Rotation 

If AL Z from (5b) were steady as well as small, then the hydrodynamic stress-free condition at the base of the 
viscous sub-layer. 




= rsin0co'(c‘) = 0, 


(6a) 



seems a fair substitute for (3d) because it filters out viscous exchange of angular momentum between the mantle and 
the main-stream. Equations (6a) and (3c) together give a free-slip main-stream boundary condition. 

Of the infinity of profiles that satisfy (6a), (3c) and (3e), consider the simple quadratic form 

2 

to(d < r < c~) = co(c~) - £~) Eco(c-) - 00 (d)] . (6b) 


Vastly more complicated forms for co(r) with much more shear and curvature tend to increase viscous dissipation 
and diffuse away more rapidly. By the well-known omega-effect, they also induce magnetic fields with much more 
shear and curvature which, in turn, tend to increase Ohmic dissipation and diffuse away more rapidly. Such 
complicated, seemingly transient, forms for co(r) are not selected because the excess dissipation in the body of the 
outer core interferes with the task at hand: cooling the core over geologic time by smaller scale thermo- 
compositional convective motions. (This is formalized in section 3). 


For uniform density throughout the core, substitution of (6b) into (5b) and setting the resulting integral for 
AL Z to zero yields, after some algebra, 


co(d) = -co(c)[- 


1 - e 5 - Q 
6 5 + Q* 


* 


(7a) 


where e = (d/c) = 0.3499 as before and 


Q* = (1 - e) 2 )[1/21 - e 5 + 5e 6 /3 -5e 7 /7) = 0.1064 
The eastward drift of the inner core from (7a, b) is 


(7b) 


G)(r < d + ) = -7.96 co(c') = 0.92 °/yr (7c) 

where the last step sets co(c') to co B (c‘) above. In this example, the solid inner core reduces the jump in to across the 
outer core; indeed, if e were zero, then co(0) would be -20co(c'). 

When uniform core density is replaced with bilinear density variation across the outer core and across the 
inner core, equation (6b) and the assumption of zero AL Z yield a modified version of (7a-c) that reduces to 

co(r < d + ) = -7.383 co(c’) = 0.85 °/yr (7d) 

for the boundary densities of Rennet et al. (1995). Increased density at depth decreases the eastward flow required 
to nun the axial core angular momentum anomaly, by about 10% in this example. 



Differential rotation profile (6b) is not only arbitrary, it does not quite satisfy (30. The shear at the inner 
boundary proves quite small, but the accompanying dissipation would not seem to aid inner core solidification and 
core cooling. In contrast, the concentration of shear induced dissipation into a thin upper boundary layer may aid 
core cooling by steepening the temperature gradient across the core-mantle interface; this speeds conductive heat 
transport into, and perhaps helps drive convective heat transport across, the deep mantle. Core motions that pump 
magnetic energy into the field outside the core and drive Ohmic dissipation in the mantle also export energy from an 

overheated core, presumably a core that cools, contracts and liberates gravitational energy in accord with 
Hamilton’s, if not Fermat’s, principle. 

3. Variational Method 

If a particular core angular momentum anomaly is stationary against perturbations in its underlying differential 
rotation, then the planet might settle into small oscillations about this state. Such a state could be energetically 
accessible if it does not cause too much dissipation D of magnetic and kinetic energy. One might thus seek a state 
wherein (i) AL is stationary against perturbations in co; (ii) there is no angular momentum exchange between the 
outer core and its bounding solids; and (iii) the dissipation in the core is stationary and minimum. This non- 
magnetic, uniformly rotating state has evidently not yet been attained, perhaps due to thermo-compositional 
convective core cooling. A non-zero differential rotation excited and maintained by smaller scale convective 
motions ought not interfere with such cooling by excesssive dissipation in the body of the core; however, a 
concentration of dissipation by the differential rotation into a thin core-mantle boundary layer and in the overlying 
mantle may aide core cooling. 

The importance attached to conditions (i), (ii), and (iii) above is thus reversed. The desired differential 
rotation profile has stationary dissipation, preferably concentrated near the outer boundary via conditions (3c, d). 
Angular momentum exchange is discouraged via inner-outer core locking conditions (3e, f); the resistive mantle is 
left unlocked, not only to increase dissipation near the core-mantle boundary, but because core viscous torque on the 
mantle turns out to be feeble and can be compensated by magnetic or topographic, if not gravitational, torque (see, 
e.g., Voorhies, 1991). With viscous boundary layers resolved, condition (6a) is not needed. 

Stationarity of D(to) + AAL z (co) is written 


5 {D + XAL Z } = 0, 


( 8 ) 



where 8 denotes perturbations in differential rotation profile and X is a Lagrange multipier. Only magnetic and 
viscous dissipation directly attributed to spherical vortex (3a) alone are included here. The stationary dissipation 
solution described below sets X = 0; however, some interesting side cases use non-zero A,. 

3.1 Magnetic and Viscous Dissipation 

For a Newtonian fluid with rate of strain e^, molecular kinematic shear viscosity v, and viscous stress py 

equal to 2pve ij - + 2pv8 i je kk /3, the viscous dissipation for solenoidal flow (3a) alone is 

D v (co) = J/J 2pv(e r 2 + e( 2)dV 
CV v v 

= J/J pv(8 r (o) 2 r 2 sin 2 0 dV 
CV 

c 

= (8tt/3) J pvr 4 (co') 2 dr. (9a) 

(see, e.g., Chandrasekahr, 1961). For a rigidly rotating inner core, the lower limit of integration in (9a) can be reset 
to d. Viscous dissipation in a high viscosity inner core can remain much less than in the outer core, provided the 
magnitude of viscous stress is similar in both regions. 

Ohmic dissipation of Amperian current (density J) in the liquid metal core of electric conductivity a 


D a = / <J _1 J 2 dV 


(9b) 


is thought to dominate viscous dissipation. By Ohm’s law for conductors moving in arbitrary magnetic field B, 


J = a(E + v x B) , 


(10a) 


an electromotive force vxB equal and opposite to electric field E causes no current and no Ohmic dissipation. 
Motions inferred in the frozen-flux approximation, wherein E + vxB is set to zero to ensure finite J despite large a, 
therefore cause no Ohmic dissipation. 

We choose the solenoidal vector potential A appropriate to the Coulomb gauge (B = VxA, E = -Vd> - d { A). 
To isolate current due to spherical vortex (3a) alone, denoted J(co), the portion of V(co)xB equal to 3 t A(co) plus 
electrostatic VO(co) must be removed. Ohm’s law for the spherical vortex is thus written 


(10b) 


J r (co) = -a[corsin0B0 + d r d>(co) + a t A r (co)] 
J 0 (to) = a[corsin0B r - r^d^oo) - a t A 0 (co)] . 


(10c) 



Faraday induction (3^B = Vxd jA) due to frozen-flux geomagnetic drift cog alone vanishes in a magnetic reference 
frame rotating at sidereal angular velocity + o)g. In such a frame, the relative angular velocity of the fluid is co^ 
by (3b); therefore, the 9 t A(co) term in (10b, c) is eliminated by replacing co with co^ 

J r (o)) = -otto^rsinOBQ + 8 r 0(co)] (iod) 

J e (co) = alco^rsinOBj. - r'^Ofco)] . (ife) 

The electrostatic term in (lOd, e) must also differentially rotate out, but with a conceivably different Ohmically 
non-dissipative vortex denoted co*(r). Clearly J(co*) must be zero, as must rJ e (co*) - 9^(0*). For uniform o, 
the latter condition applied to (lOd, e) gives 

o[9 r (cD*r 2 sin0B r ) + dQ(co*rsin0B0)] = 0, (Ha) 

or, for solenoidal B, 


o[(co*) ' AinOBj. - = 0. (Ub) 

For to to be non-zero, shearing of the field must balance advection of the non-axisymetric field. 

A magnetically dissipative vortex co^ - co* violates (11), generates non-zero meridional currents (10), and 
thus induces azimuthal field from meridional field (the famous omega-effect). Substitution of (lOd, e), after 
removing non-dissipative co*, into (9b) yields 


D ct (co) = J / / 0(01^ - co*)Vsin^0(B r ^ + Bq^) dV . 


The definition of weighted mean meridional magnetic pressure. 
In n 

P m ( r ) s (3/87C) / J (2p.)'l(B r 2 + B0^)sin^0d0d(j) , 
0 0 


and (3b) allow (12a) to be rewritten as 


c 

D a (co) = (8tt/3)J 2P m ir 1 (co - co B - co*) 2 r 4 dr, 
0 


(12a) 


(12b) 


(12c) 


where q is magnetic diffusivity Ota)’ 1 . If to equals co B + co* in a magnetically locked inner core, the lower limit of 
integration in (12c) can be reset to d. Elsewhere, such balance minimizes Ohmic dissipation due to the vortex alone. 
Note (12c) excludes the inner product of other currents within the core, notably sources of P m , with J(co). 


12 . 


When co on the left of (1 lb) is replaced with magnetically dissipative co^ - co*, the resulting term is not 
generally zero. A non-zero integral of the square of this term is thus a sure sign of, if not directly proportional to, 
Ohmic dissipation. The imbalance between shear and advection of the magnetic field by a dissipative vortex 
suggests a possible alternative expression for Ohmic dissipation 

D a (w) ~ JJ J cKco^ - co*) 'r 2 sin0B r - (co^ - ©Vd^B^dV . (13) 

If the advective term proportional to were omitted, then (13) would resemble (9a) with an ’effective viscosity 
p%’ equal to crr^Bj. 2 . With a ~3xl0~* S/m and B f 2 » (3.3 gauss) 2 averaged over I r| =c, px is about 4x10^ Pa 
s and x is about 4x 10 2 m 2 /s. The latter value is about 14 orders of magnitude larger than the molecular kinematic 
shear viscosity of the liquid metal. This may help explain differences, highlighted by Lumb & Aldridge (1991), 
between ’effective viscosities’ obtained presuming a non-magnetic core and the viscosity of the liquid metal. 
Appendix A shows how such oversimplification of D a (co) destroys boundary layer structure yet, with non-zero X, 
leads to seemingly reasonable values for inner core eastward drift. 

3.2 Variational Calculus and Illustrative Solution 

With AL Z from (5b) and D(co, co') equal to the sum of magnetic dissipation (12c) and viscous dissipation (9a), 
the stationarity condition (8) is written 
c 

5 { (87C/3) / [2P m rf^co - cog - co*) 2 + pv(co') 2 + ?ipco]r 4 dr} = 0. (14) 

d 

The lower limit of integration has been reset to d for a rigidly rotating, magnetically locked inner core. Stationarity 
requires the integrand of (14), denoted f(co, co ; r), to satisfy Euler’s equation, which sets 3f/dco equal to (df/dco ' ) ' 
(see, e.g., Marion, 1970) and reduces to 

co'' + [ln(pvr 4 )] ' co ' - 2P m (pvr|)' 1 (co - co B - co*) = X/2v ( 15 ) 

To simplify (15), omit variations in p, v, and T) with radius and set X to zero. For such stationary dissipation. 


co'' + (4/r)co ' - 2P m (pvTi) _1 (co - co B - co*) - 0 . (16) 

Equation (16) also describes a balance between azimuthal Lorentz and viscous forces (see Appendix B). The 

general solution to (16) is the sum of particular solution C 0 p and homogeneous solution Particular solution C 0 p = 
* 

cog + co satisfies 

co p " + (4/r)co p ' = 0 


/ / 


(17a) 



and is 


Wp(r) ~ co 0 + B(c/r) 3 , (17b) 

where co 0 and B are integration constants. Primary spherical Couette flow (17b) also describes a non-magnetic case 
(Cartwright et al„ 1996; B. Bills, 1998); non-zero X adds a quadratic to co p (r). Homogeneous solution % satisfies 
5" + (4/r)£' - 2P m (pvr|)- 1 £ = 0. (18a) 

The functional form of P m (r) is needed to solve (18a). The detailed form does not matter much for the problem at 
hand, provided P m (r) does not vary by more than eight orders of magnitude across the outer core, because 
2P m r 2 (pvr|)- 1 is enormous and £ is therefore small except in thin boundary layers. 

If the meridional magnetic field originates in the core, then positive P m (r) should tend to increases with depth 
below c, perhaps reaching a maximum and then falling to modest values. The illustrative case simply takes P m (r) to 
be Kc 2 /! -2 , so solutions to (18a) are of power law form 

% = C(r/c)P + D(r/d)l (i 8b ) 

where 

p = (1/2) {-3 + [9 + (SKc^pvn)] 1 / 2 } 

» (2Kc 2 /pvr|) 1 ^ = 2xl0 7 . (Igc) 

q = (1/2) {-3 - [9 - (8Kc 2 /pvr|)]^} - -p 

The value of p in (18c) follows from; a 5 gauss meridional field at radius c = 3.48 Mm, hence K~ 0.1 Pa; p~ 10 4 
kg/m , the molecular kinematic shear viscosity v of the molten metal, which is well-known to be about 3xl0‘ 7 m 2 /s 
(Poirier, 1988); and o « 3x 10 5 S/m as above. Because p is so large, the first and second terms in (18b) respectively 
dominate the core-mantle (outer) and inner-outer core (inner) boundary layers. The outer boundary layer scale 
height I £/£ J is 8 = c/p = 20 cm (about 2.5 times the Ekman depth). The inner boundary layer scale height is e8. 
If p m ( r )EP VT l]" 1 were decreased by eight orders of magnitude, 8 would increase to 2 km, which remains thin 
compared with 2.2 Mm of molten iron. The illustrative case is therefore robust. 

The general solution of (16) is the sum of (17b) and (18b) 


co(r) = co 0 + B(c/r) 3 + C(r/c)P + D(r/d)d . ( 19 ) 

Spherical vortex (19) exerts no net azimuthal force density in the outer core in so far as viscous diffusion of its 
momentum balances Lorentz drag. 



3.3 Boundary Conditions for Magnetically Locked Inner Core 

The four parameters (co 0 , B, C, D) in (19) are precisely those needed to match bounday conditions (3c-f). The 
contribution to ©(c) from D(c/d)^ = D(c/d)**q ** De**(2xl0 7 ) proves negligibly small, so the no-slip core-mantle 
interface condition (3d) on (19) implies 


co. 


B - C. 


J o ~ ----- ( 20 a) 

Westward drift condition (3c) is applied at c' ~ c - 105. To an accuracy of [(c - 105)/c]P~ 10"^, this gives 


©(c“) = co Q + B = - C . 

Similarly neglecting the small contribution from CeP to co(d), 

co(d) = B[(c/d) 3 - 1] - C + D . 

The top of the inner boundary layer is taken to be d + = d + lOeS, so 


(20b) 


(20c) 


o)(d + ) - B[(c/d) 3 - 1] - C . (20d ) 

By subtracting (20c) from (20d), we see that the jump in © across the inner boundary layer is -D. Thus I D i is small 
when © ' (d) is zero, as required by condition (3f). The latter is intended to reduce dissipation in the inner-outer core 
transition region that does not immediately aid core cooling. Because CpeP « C, this condition gives 


CO '(d) = - (3B/c)(c/d) 4 + (Cp/c)(d/c)P _1 + (Dq/d) 

= (Dq/d) - 3B(c/d) 4 = 0 , 

hence 

D = 3q _1 B(c/d) 3 . 

Recalling q = -2xl0 7 , (20f) confirms I D I « I B I . Substitution of (20a, b, f) into (19) gives 


(20e) 


(20f) 


w(r) = ©(c")[l - (r/c)P] + B[(c/r) 3 - 1 + (3/ q~) (c/d) 3 (r/d) °1] ( 21 ) 

To relate B to ©(c - ), and thus assign a numerical value to co(d) via (3c) and (20d), it is assumed that AL Z in 
(5b) is zero. (Alternatively, AL Z might be calculated from a value for ©(d) or constrained via non-zero X). With p 
already treated as uniform and omitting the tiny contribution to AL Z from the boundary layers, this requires 
c' 

0 = ®(d)d 5 /5 +J {©(c - ) + B[(c/r) 3 - 1] }r 4 dr . 
d + 


(22a) 



(22b) 


The elementary integration, with 00 (d) being co(c') + B[(c/d) 3 - 1) from (20b, d), gives 
B = -2co(c')[3 -e 2 ]" 1 . 

Substitution of (22b) into (21) yields the solution 


co(r) = co(c-){ [l-(r/c)P] + 2(3 -e 2 )' 1 ^) 3 - 1 + (3/q)(c/d) 3 (r/d)q] } . (22c) 

Clearly (22c) depends wholly upon oo(c'). If co(c') were presumed to be zero, then co(r) would vanish; if | co(c") I 
were presumed to vastly exceed I oo B ( C -) I , then co(r) would be quite large. Geomagnetic inversions indicate an 
intermediate case wherein co(c’) is a large fraction of w B (c'). 


By (22c), or by combining (22b) with (20d) and (20b), 

_ , 2(e 3 - 1) 

co(d) = co(c‘)[l + -1 

3e 3 -e5 

= -co(c')[ 13.31] . 


(23a) 

(23b) 


In the diffusive extreme, co(c") and co(d) would be zero. In the frozen-flux extreme, co(c") would be co B and, using 
cog from (3b), 


co(d) = -w b (c-)[ 13.31] = 1.54°/yr . (23c) 

Equation (23c) offers a soft upper bound on inner core eastward drift. 

The shear across the core-mantle boundary layer in the illustrative example is 

co'(c) = q + C/c - 2B/c « -(0(0/5 . (24) 

The feeble viscous torque exerted by (24) on the mantle must be compensated by magnetic torques on the resistive 
mantle or magnetically permeable crust; however, even with the frozen-flux value for co(c"), a very weak toroidal 
magnetic field at the core-mantle interface, albeit of magnitude no less than 8 nT, could suffice (Voorhies, 1991). If 
it were not compensated, the resulting secular increase in the length of the day of roughly 5.5 |iis/decade, or less than 
one hour per 4.5 Gyr, would seem far smaller than other effects omitted above. 

Both viscous and magnetic dissipation due soley to spherical vortex (22c) are concentrated in the core-mantle 
boundary layer and sum to a tiny fraction of the geothermal flux. Viscous dissipation (9a) in the core-mantle 
boundary layer due vortex (22c) alone would be about (47cc 4 pv/35)[w(c')] 2 , or 40 kW. Magnetic dissipation (12c) 
due to this vortex alone is due solely to the homogeneous solution (18b). For the form of P m adopted above, the 



portion of this disspation occuring in the core-boundarly layer would be about (47rc 4 Kpa/3S)[co(c-)] 2 , or about 0.5 
MW. These values seem very small; however, if other motions of the outer core dissipated magnetic and kinetic 
energy at the same rate per unit volume, then total disspation in the outer core would be about 3 TW. Dissipation of 
this magnitude might arise if the other motions continually entrain and mix thin boundary layers into the body of the 
outer core or, perhaps more simply, have a very small scale in one direction. 

4. Summary 

The idea that geomagnetic westward drift indicates a slight convective leveling of the momentum gradient 
within Earth’s core was pursued in search for a differentially rotating mean state, upon which various oscillations 
and secular effects might be superimposed. The variational calculus of stationary dissipation applied to a spherical 
vortex within Earth’s core leads to an inhomogeneous second order differential equation for the differential rotation 
profile w(r). For a magnetic, liquid metallic outer core (P m r 2 /pvT| » 1), the bulk of the vortex is separated from 
the bounding solids by 20 cm thin magneto-viscous boundary layers. 

Because part of the motion may induce no electric current, hence no Ohmic dissipation, four boundary 
conditions are required instead of two. The four conditions imposed on this form are: (i) no-slip at the core-mantle 
interface, (ii) geomagnetically estimated bulk westward flow of up to 0.12°/yr at the base of the thin core-mantle 
sub-layer, (iii) no-slip at the inner-outer core interface, and (iv) the hydrodynamically stress-free inner core boundary 
condition appropriate to a magnetically locked inner core. To compute the eastward drift of such an inner core from 
geomagnetic westward drift, it is assumed that the axial core angular momentum anomaly is zero. 

The eastward drift of the inner core resulting from this analytic exercise is at most 1.5°/yr. This value is 
within the 0-3°/yr range predicted by some numerical simulations and inferred from seismological data. It is 
remarkable that a simple, kinematic spherical vortex model featuring liquid metallic viscosity should show any 
similarity to either. Values temporarily exceeding 1.5°/yr might be explicable in terms of decadal oscillations. 
Provided the axial angular momentum anomaly of the core is not too large, there appears to be no serious conflict 
between geomagnetically inferred westward drift at the top of the core, eastward drift at depth, and seismologically 
inferred eastward drift of the inner core. The general agreement by no means reduces the importance of inner-outer 
core-mantle oscillations, differential precession, differential tidal despin, or differential spin up due to inner core 


solidification. 
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APPENDIX A 

If d<j) B <J> were omitted and heterogeneous ar 2 B r 2 parameterized as a laterally homogeneous ’effective 
viscosity’ p(r)x(r), then (13) would reduce to 

c 

D a » (87i/3) / px r 4 (co ' ) 2 dr . ^ 

The similarity between (Al) and (9a) provokes comment on dynamic shear viscosity pv as a parameterization 
of momentum transport by electromagnetic molecular interactions ("collisions'’). Hydrodynamic stress (p~ = -pSjy + 
p ij) ma y be vie wed as a macroscopic average over intermolecular electromagnetic (Maxwell) stress <T i j>. The 
reference stress T°jj is due to a tangle of fluctuating intermolecular fields (e, b) originating in the molecules. A 
macroscopic average (e.g., over a cubic millimeter and a milli-second) amounts to an isotropic cohesive equilibrium 
<T ij 0> = -p5ij. Applied macroscopic strain rate to'rsinG shears the intermolecular field within a liquid, arguably 
inducing a contribution to b^ equal to co ' rsin0Kb r , where k depends upon details of the interactions. The 
contribution to <T r ^> is <b r b ( j ) >/2p 0 , or (rsin0/2p o )<Kb r 2 >co ' . In the absence of macroscopic fields, <b> remains 
zero; moreover, terms analogous to those omitted in going from (13) to (Al), specifically <0 <J) b ( j ) ) 2 > and <b r dfy>, 
would seem as negligible as <5 ( j,b ( | ) > when the applied strain rate is uniform over the averaging 4-volume. If the 
contribution to <T r( j ) > is identified with Py, then pv might be identified with <Kb r 2 >/4ju 0 . The fact remains that 
interactions establishing viscosity involve intermolecular collisions, while macroscopic B is due to Amperian 
currents regulated by conduction electron collisions; therefore, k differs from a and v differs from %. Large values 

of core effective viscosity’ have little or nothing to do with viscosity, but may have much to do with Loren tz forces 
and Ohmic dissipation. 

If D(co) were treated as the sum of (Al) and (9a), then (8) would amount to 
c 

0*6 (8jt/ 3) J ^ [p%(co ' ) 2 + Ap 0 (r)co(r)]r 4 dr ^ 2 ) 


because v « %■ Stationary would require the integrand of (A2) to satisfy Euler’s equation. 



to' ' + Un(p%r 4 )] ' co ' - X/2% = 0. ( A3 ) 

The first and second integrals of (A3) are straightforward when p(r) and %(r) are, respectively, approximated by 
uniform constants. The result is written 


co(r) = 0 )! - (C 1 /3px)r“ 3 + (A/2 O^r 2 (A4) 

where Cj and coj are the two constants of the two integrations. In (A4), primary spherical Couette flow is the 
homogeneous solution and the parabolic term on the far right is the particular solution for non-zero A. The constrast 
with (19) is remarkable; moreover, (A4) shows no sign of the boundary layers and, even with non-zero A, can sate 
but three of four boundary conditions. 

Stress Free CMB vs. ICB: If co (c') were zero, then (A4) would imply 


(A/X) = -KXCj/p^c- 5 . 

When (A4, A5a) are used in (5b) and AL Z is set to zero, again for uniform density, one obtains 
(Cj/PX) = C0jC 3 /5Q , 

where Q depends only on e and is about 0.2259. It has also been shown that 


(A5a) 


(A5b) 


coj - co(c’) 6Q/(6Q - 1) . (A5c] 

Subsitution of (A5a-c) into (A4) and evaluation of the result at inner radius d yields 

co(d) = co(c') [6Q/(6Q - 1)] { 1 - (5Q)' 1 [(3e 3 )' 1 + e 2 /2]} 

= - C0(O (22.68) = 2.62°/yr , (A5d) 

where the last step uses the frozen-flux value for co(c") from (3c). One also obtains an coj of -2.44xl0" 10 rad s" 1 
and C/px = -9.1xl0 9 m 3 /s. 

The foregoing oversimplification of D a (co) leaves some shear at the inner core boundary; indeed, from (A4 - 
A5d) it follows that 


co ' (d) = co(c) c 4 (6/5)(6Q - 1)" 1 Pe" 4 - e] 

= 224.2 c _1 co(c') = -28.3 co(d)/d . (A6) 

Such shear could curl field lines threading the inner core and generate non-irrotational currents, magnetic diffusion. 
Ohmic dissipation that might inhibit inner core solidification, and restoring Lorentz torques. If it were not 



compensated by other torques, the viscous stress would tend to spin down the inner core and spin up the outer core, 
albeit extremely slowly. 

If the hydrodynamic stress-free condition were instead applied at the inner boundary, then a slower inner core 
eastward drift of -1.847co(c), or about 0.21°/yr, would result. The non-zero shear on the main-stream below the 
core-mantle boundary, to (c ), would be about be about 7.08 co(c')/c. If uncompensated, this would tend to spin up 
the core and spin down the mantle. 

APPENDIX B: Force Balances 

The azimuthal Lorentz force density implied by (lOa-e) is J r (a))B 0 - J 0 (co)B r , is -a(co - cog ci)*)rsin9(B r 2 + 
B 0 2 ), and is equivalent to magnetic friction -Cf(V^ - Vg - V*) with drag coefficient Cyr = a(B f 2 + B 0 2 ). The 
azimuthal viscous force density from the divergence of p r ^ is, assuming uniform pv, pvsin0(rco ' ' + 4co ‘ ), or 
pvr (r co ) . There is no azimuthal Coriolis force on motions v^; pressure forces are scaloidal; and buoyancy 
forces are mainly radial. A balance between Lorentz drag and viscous diffusion would have 

pvsin0(rco "+4(0') - a(co - cog - <o*)rsm0<B r 2 + b 0 2 > _ 0 ^ 

Only the azimuthal mean value of (B f 2 + B 0 2 ), denoted by angle brackets, appears in (B 1) because other currents J - 
J(co), and other motions v - V(co), in the core are excluded from the analysis of the vortex alone. When (Bl) is 
multiplied by sin 0d0d<{> and integrated over a sphere, the result is identical to (16). Spherical vortex (19) is thus 
force-free in the limited sense that the viscous diffusion of its momentum is balanced by the Lorentz force caused 
by the current it generates from the ambient field. 
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